Charge Current Density from the Scattering Matrix 
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A method to derive the charge current density and its quantum mechanical correlation from 
the scattering matrix is discussed for quantum scattering systems described by a time-dependent 
Hamiltonian operator. The current density and charge density are expressed with the help of 
functional derivatives with respect to the vector potential and the electric potential. A condition 
imposed by the requirement that these local quantities are gauge invariant is considered. Our 
formulas lead to a direct relation between the local density of states and the total current density 
at a given energy. To illustrate the results we consider, as an example, a chiral ladder model. 



( Introduction } The scattering matrix gives an impor- 
tant starting point in descriptions of quantum transport 
phenomena. It connects the incoming current amplitudes 
to the outgoing current amplitudes, and is calculated 
from the Hamiltonian operator by the M0ller operator 
or the Green function method Q). We can obtain infor- 
mation about global characteristics in systems from the 
scattering matrix directly. For example, the Landauer 
formula gives a method to calculate the conductance from 
transmission amplitudes as components of the scattering 
matrix The Friedel sum rule connects the density 

of states to the scattering matrix |fi§|, so that it is in 
principle possible to calculate any equilibrium statistical 
mechanical quantity from the scattering matrix |10|. We 
can also derive an expression of the persistent current in 
open conductors caused by a magnetic flux p| , or the 
shot noise [12jlq] from the scattering matrix. 

It shouloibe emphasized that the scattering matrix 
even gives information about local characteristics in sys- 
tems. We consider one-particle and time-dependent 
Hamiltonian systems. First, the charge current density 
J n (x,t) at position x and time t generated by the par- 
ticle incident in a state with the quantum number n is 
connected to the scattering matrix S — (S nn i ) as 



(1) 



where A(x,t) is the vector potential and c is the velocity 
of light. In this Letter we call this quantity the "injected 
current density" . Eq. (Q) is the main result of this Let- 
ter. Second, the probability density p n (x,t) of the state 
caused by the incident particle with the quantum number 
n is given by 



pn(x,t) 



2ni ^ " " 



6U(x,t) 



(2) 



with the potential U(x,t). The probability density 
p n (x, t) also has been called the "injectivity" and can be 
regarded as the density of states with a preselection of 
the incident channel. Eq. (||) in the time-independent 
Hamiltonian case has been used in treatments of the 



electron-electron interaction using a self-consistent po- 
tential [p^-p7[. In this Letter we give a new derivation 
of this formula and generalize it to the time-dependent 
Hamiltonian case. 

Our technique to derive Eqs. ([!]) and (j^) also can be 
used to arrive at quantum mechanical correlation func- 
tions of local quantities from the scattering matrix. For 
example we show that the quantum mechanical correla- 
tion C^ u) {x,x';t) between the /i-th charge current den- 
sity component at position x and the v-th. charge current 
density component at position x' at the same time t is 
given by 



ss„ 



5A^{x,t) SA^(x',t) 



(3) 



where A M (x,t) is the /i-th component of the vector po- 
tential. 

Our formulas include the vector potential and the po- 
tential explicitly, so we have to discuss the gauge invari- 
ance. We show Eqs. (Q - 1|) and (0) to be gauge invariant, 
and derive conditions which should be satisfied by locally 
gauge invariant quantities. 

In the time-dependent Hamiltonian system an energy 
of an outgoing particle can be different from the energy 
of the corresponding incoming particle, and the scatter- 
ing matrix element S nn > describes a transition to such 
a different energy. On the other hand, if we consider 
only time-independent Hamiltonian cases, the problem 
becomes simpler, because the scattering matrix is de- 
composed into the scattering matrices restricted to the 
energy shells. Using this feature we obtain a relation be- 
tween the local density of states and the total current 
density defined by the sum of the injected current den- 
sity J n (x,t) with respect to the suffix n satisfying the 
condition E n = E at energy E. 

As a simple example we investigate a ladder model 
with a directionality, namely a chiral ladder model, 
termed a "quantum rail road" in Ref. [Q. We verify 
the formula ([!]) for this model calculating separately the 
current density and the scattering matrix. 
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( Injected current density and injectivity ) We start 
from the time-dependent Hamiltonian operator H(t) in 
the quantum scattering system. The dynamics of the 
system is described by the Schrodinger equation using 
this Hamiltonian operator. We decompose the total 
Hamiltonian operator H(t) into the asymptotic Hamil- 
tonian operator Hq and the scattering operator -ffi(i); 
H(t) = Hq + Hi(t). The operator Hq is chosen to be 
the Hamiltonian operator which describes incoming and 
outgoing particles in the asymptotic regions. We as- 
sume that Hq is a time-independent operator determined 
uniquely. Besides, the system described by the Hamilto- 
nian operator H(t) or Hq is assumed to have no bound 
state. Below we use the coordinate representation for any 
operator and take x as the coordinate of particles. Un- 
der these conditions the scattering matrix element S nn i 
is given by 

S nn , = lim lim dx$ n {x)*U(t 2 ,ti)$ n >(x) (4) 

t-2 — > + 0O t\ — > — OO / 

where U(t2,ti) is the time evolution operator 
e^p{iH t 2 /h}-fe^p{-i J^dtH(t)/h} -exp{-i Hoh/h} 

in the interaction picture with T being the positive time- 
ordering operator, and <&„(a;) is the eigenstate of the 
operator Hq corresponding to the energy eigenvalue E n . 
Here the limits t\ — > — oo and t 2 — > +00 are defined 
by \im t ^ ±00 X(t) = lim^ +0 (±e) J Q ±oc dt'e^' X{t') for 
any function X(t) of t 0. The set {$ n (x)} n of the 
eigenstates of the operator Hq is chosen to satisfy the 
orthonormality condition and the completeness relation. 
The scattering matrix S = (S nn i) is shown to be an 
unitary matrix; S^S = SS^ 1 = I. 

For simplicity wc consider the one-particle system. 
The total Hamiltonian operator H (t) and the operator 
Hq are represented as (—ihd/dx — qA(x,t)/c) 2 /(2m) + 
U(x,t) + U Q {x) and -(h 2 / (2m))d 2 / dx 2 + U (x), respec- 
tively, where m is the mass of the particle, q is the charge 
of the particle, A(x, t) is the vector potential, and U{x, t) 
is the external potential (plus the induced potential by 
the interaction of particles), and Uq(x) is the confinement 
potential. We consider a functional derivative 5 {t> S nn ' of 
the scattering matrix element with respect to the vector 
potential A(x,t) or the potential U(x,t) at time t, us- 
ing the notation 5 {n = S/SA(x,t) or S/5U(x,t). Using 
the expression (^) of the scattering matrix elements we 
obtain a general relation 

1 / r°° \ (t) 

^-Y, S n>n(S it) Sn' n ) = -( dt' 6™ H^t') ) .(5) 

n t \ J — oo I n 

Here the notation ( • • • )< l t) means the expectation value 
taken with the scattering state of the quantum num- 
ber n, namely (X)%> = J dx^ n (x, t)*X^ n (x, t) 
/(2irh) for any operator X, where ty n (x,t) is a 
solution of the Schrodinger equation using the to- 
tal Hamiltonian operator H(t) and is defined by 



¥„(as,t) = lim tl __ M fexp{-ij; t l dt / H(t')/h}- 
exp{—iE n ti/h}$ n (x). The derivation of Eq. (||) is 
given by using the completeness relation of the set 
{§ n (x)} n , the property U(h,ti) = U(t 2 ,t)U(t,ti) of the 
time evolution operator and the relation 5 {t) U(t 2 , h) = 
U(t 2 ,t) ■ exp{iH t/h} ■ {{-i/K) f^dt'S^H^t')} ■ 

cxp{-iH t/h} ■ U(t,ti) in ti < t < t 2 . Eq. |j) is a 
key result of this Letter. We introduce the injected cur- 
rent density J n (x,t) and the injectivity p n (x,t) as 

J n ( X ,t) EE ( J(x,t) }™ , p n (x,t) EE { P (X) )W (6) 

using the charge current density operator J(r,t) = 
qv(t) ★ 8{x — r) with v(t) being the velocity oper- 
ator (—ihd/dx — qA(x,t)/c)/m and the multiplica- 
tion * being the symmetrized product, and using the 
probability density operator p(r) = 5(x — r). Eqs. 
(0) and ^ are derived from Eq. (jEJ), using the 
relations J^° oo dt'6Hi(t')/5A(x,t) — —J(x,t)/c and 
dt'6H 1 (7)/5U{x, t) = p(x). 

( Current density correlation ) We consider functional 
derivatives 5 { p S nn >, j — 1, 2 of the scattering matrix ele- 
ment with respect to the vector potential A(x,t) or the 
potential U(x,t) at time t. Using Eq. (||) we obtain 
another general relation 

n' 

= ( dt'S^H^ (J™ dt"8fH,{t")^ Y . (7) 

Using Eq. (^) we obtain Eq. (p|). Here the quantum 
mechanical current density correlation C^*" 5 (x, x'; t) is 
defined by 

C<r>(x,x';t) ee ( P»(x,t)*jM(x',t) ) m (8) 

\ In 

where J l,i) (x,t) is the p-th component of the charge cur- 
rent density operator. It should be noted that the aver- 
age ( • ■ • }^ ) taken in the correlation (@) includes only the 
quantum mechanical average, but does not include the 
thcrmo-statistical average. 

In similar ways we can obtain other quantum me- 
chanical correlation functions ( J (v> (x, t) * p(x') )[ i t) and 
( p(x) * p(x') )[ l t> from the scattering matrix. 

( Gauge invariance ) The gauge transformation of 
the electric potential (j>(x,t) and the vector potential 
A(x,t) is represented as (j>(x,t) — > 4>(x,t) + A(j>(x,t) 
and A(x,t) -> A(x,t) + AA(x,t) with Acp(x,t) = 
— (l/c)dip(x,t)/dt and AA(x,t) = d(p(x,t)/dx using a 
function ip(x,t) of x. The electric and magnetic fields 
are invariant under the gauge transformation. Using that 
the Schrodinger equation is gauge invariant, we can show 
that the p-th component J!^ ] (x, t) of the injected current 
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density is also gauge invariant, meaning that the formula 
(|l|) is gauge invariant. 

The fact that the injected current density is gauge in- 
variant is also represented as 



dt' 



,i x ' { yfffoff a^x',^) 



(x',t>) 
SA(x',t') 



AA(x',t')} = 0. (9) 



A change 5<f>(x, t) of the electric potential <f>(x, t) modifies 
the potential U(x,t) as 6U(x,t) = q5<fi(x,t). Moreover 
Eq. (^) should be satisfied for any function ip(x,t) which 
is zero at the boundary of the integral of the left-hand 
side of Eq. (^) . Using these facts and the partial integrals 
in Eq. (ph we obtain 



q 8 8J}?>(x,t) 
c dt' 8U(x', t>) 



d 

dx' 



5A(x',t') 



= 0. 



(10) 



Eq. (10) represents a condition for the injected current 
density imposed by its gauge invariance. Similarly, the 
formulas (0) and (0) are gauge invariant, and the injec- 
tivity p n (x,t) satisfies an equation similar to Eq. (|l0|). 

( Total current density and local density of states in 
the time-independent Hamiltonian system ) Now, we 
consider the time-independent Hamiltonian case. In this 
case the injected current density and the injectivity are 
independent of time t. Moreover the state ^ n (x,t) used 
to define the average ( • • • )[ l t) can be replaced by the 
state *„(cc) = {1 + G(E n )Hi}<& n (x) using the Green 
operator G(E) = lim e ^ +0 (-E — H + ie) -1 . The state 
^ n (x) is the eigenstate of the total Hamiltonian oper- 
ator H corresponding to the eigenvalue E n and satisfies 
the Lippmann-Schwinger equation. One of the important 
features in the time-independent Hamiltonian system is 
that the scattering matrix element S nn > takes a non-zero 
value only in the case of E n — E n i . Therefore, in the 
formulas ([!]) and (||) we can replace the full scattering 
matrix element S n > n with the matrix element Si>i(E n ) 
of the scattering matrix S(E n ) restricted to the energy 
shell of energy E n , and exchange the sum of the quantum 
number n' with the sum of the channel number I'. 

It is important to note a difference between the 
time-dependent and time-independent Hamiltonian cases 
in our formulas. In the time-dependent Hamiltonian 
case the functional derivative 5S nn ' /5A(x, t) is given 
by 5S nn > = J dx J dt (SS nn > /SA(x,t))5A(x,t). On 
the other hand, in the time-independent Hamiltonian 
case the functional derivative SSw (E)/SA(x) is given 
by 6S W (E) = J dx{8S w (E)/SA(x))SA(x), which does 
not include the time integral like in the time-dependent 
Hamiltonian case. This difference appears in the cur- 
rent density correlation formula (^), where a simple ex- 
change of the scattering matrix element S nn i and the 
vector potential A(x,t) with the on-shell scattering ma- 
trix element Sw(E) and the vector potential A(x), re- 
spectively, is not allowed to obtain its time-independent 



Hamiltonian version. However we can obtain the time- 
independent Hamiltonian versions of the injected current 
density formula (Q) and the injectivity formula (^) by 
such a simple exchange. This technical point is better 
considered in a separate paper. 

We proceed to consider the total current density 
J(E, x) and the local density of states p(E, x) defined 
by 

J(E,X)=J2 Jn(x), P (E,X)=J2 Pn(x). (11) 

n n 
(e„=e) (_E„=_e) 



Eqs. (|l|) and <g) lead to 
7r dA(x) 



p(E,x) 



1 S9 f (E) 
"tt SU(x) 



(12) 



where 9f(E) is called "Friedel phase" Q and is de- 
fined by 9 f (E) = (l/(2i)) \nDet{S(E)}. It may be noted 
that using the first equation in Eq. ([l2]) the relation 
J(E,x) = —J(E,x)\a-+-a is obtained under the con- 
dition S nn > = S n > n \A^-A- The second formula in Eq. 
( |l2| ) have already been used in some works |l^- [l7| , ^o]| . It 
follows from Eq. (jl^) that 



5J{E,x') 
SU(x) 



8p(E,x) 
' SA(x') ■ 



(13) 



This equation shows a connection of the current density 
to the local density of states. 

The total current density J(E, x) takes a non-zero 
value in some important phenomena, such as edge cur- 
rents in the quantum hall effect, persistent currents 
caused by a magnetic flux and so on. Earlier work |Tl| ] 
obtained the total persistent current from the flux deriva- 
tive of the scattering matrix for a ring connected to a 
lead pH . For a purely one-dimensional ring the func- 
tional derivative with respect to the vector potential can 
simply be replaced by the derivative with respect to the 
flux. In such a special case the first formula in Eq. jl2] ) 
leads to the same expression as the earlier work. 

( Example ) As a specific example we consider a chiral 
ladder model threaded by a weak magnetic field. In this 
model particles move on the both legs of the ladder only 
in one direction (See Fig. |l]). This may be regarded as a 
model of two edge channels at one edge of a quantum Hall 
bar with impurities, by which electrons transfer from an 
edge channel to another channel. 
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FIG. 1. Chiral Ladder Model. 
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We assume that the ladder has one-dimensional legs 
each of which has one channel. We introduce the scat- 
tering matrix T (j) which connects the current amplitude 
to the left of the j-th rung to the current amplitude to the 
left of the j+l-th rung. The dependence of the scattering 
matrix T <J> on the vector potential in the legs is 



tg exp{#!f7<M t { H exp{i^V0o} 



(14) 



where <fio = Tic/q and (pf (02 ) i s the integral of A(x) 
over the upper (lower) leg between the j-th and the j+l- 
th rungs with A(x) being the vector potential element 
in the direction of the legs at position x. Here the ma- 
trix t (j) = (tni) is the corresponding scattering matrix in 
the case that the vector potential in the legs is zero. In 
this model the scattering matrix T u) is the same as the 
corresponding transfer matrix, so the scattering matrix 
S u} = (Sin) of the sub-system consisting of the first j 
number of rungs is given by T^'T' 3-1 ' • ■ • T (1) . 

We consider the injected current density J{ ] J tl) (J^y) 
in the upper (lower) leg between the j-th. and the j + l- 
th rungs, which is caused by the incident current a M 
shown in Fig. |l|. This current is represented as J$ lV ) = 
qpQv\S { J^\ 2 where v is the particle velocity in the up- 
per (or lower) leg and po is the local density of states 
l/(27r?i|w|) in the one-dimensional perfect wire. 

Now we connect the scattering matrix S (N> of the sys- 
tem consisting of N number of the rungs (N > 1) to 
the injected current density J<^,„). Noting unitarity of 
the matrices T u '\ j'=N,N-i,---,i+i and using the relation 



dij'Svv', where x 



„(i) 



(rr!, ) is a point 



in the upper (lower) leg between the j-th and the j+l-th 
rungs, we obtain (1/(2™)) £li S<™ ■ 5S™/8A(xP) = 
(1/c) Jt,] U ), which is just the injected current density for- 



mula (|l|) in the time-independent Hamiltonian case. 

( Conclusion and remarks ) In this Letter we have 
discussed formulas to derive the charge current density, 
the charge density and their quantum mechanical current 
density correlations from the scattering matrix in one- 
particle and time-dependent Hamiltonian systems. The 
gauge invariance requires Eq. ( [To|) which has to be satis- 
fied by these local quantities. Using our formulas we ob- 
tained a relation between the local density of states and 
the total current density produced by incident particles 
at an energy in time- independent Hamiltonian systems. 
Specifically we verified the current density formula for a 
chiral ladder model. 

In this Letter for simplicity we considered one-particle 
systems only. However our technique to derive Eqs. (|5|) 
and (0) can be used to obtain some generalizations to the 
formulas including many-particles' effects. For instance, 
we can generalize Eq. (pi) to the inelastic scattering cases 
by dynamical scatterers with no charge. Eqs. (ra) and 
(pi) also suggest that we can derive formulas for other 
local physical quantities from functional derivatives of 
the scattering matrix with respect to other local fields. 
For example, if the potential U(x,t) includes the term 
— qs ■ B(x,t)/2 as an interaction effect of a spin s with 



a magnetic field B(x,t), then we can calculate the local 
expectation value of a spin component from a functional 
derivative of the scattering matrix with respect to the 
magnetic field. 

As another approach to the char ge c urrent density we 
can use the linear response theory |22|23[|. The connec- 



tion of the scattering theoretical approach to the linear 
response theoretical approach in the description of local 
quantities is left as a future problem. 
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